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ABSTRACT 

We investigate the pulsational properties of Pre-Main-Sequence (PMS) stars 
by means of linear and nonlinear calculations. The equilibrium models were 
taken from models evolved from the protostellar birthline to the ZAMS for 
masses in the range 1 to 4 M . The nonlinear analysis allows us to define the 
instability strip of PMS stars in the HR diagram. These models are used to 
constrain the internal structure of young stars and to test evolutionary models. 
We compare our results with observations of the best case of a pulsating young 
star, HR 5999, and we also identify possible candidates for pulsational variability 
among known Herbig Ae/Be stars which are located within or close to the 
instability strip boundaries. 

Subject headings: stars: formation - stars: evolution - stars: pre-main- sequence 
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1. Introduction 

Young stars are characterized by a large degree of activity. Winds, jets and outflows are 
manifestations of the interaction of the stars with the circumstellar environment in which 
they are embedded. Similarly, both T Tauri and Herbig Ae/Be stars show photometric and 
spectroscopic variability on time scales of minutes to years, indicating that photospheric 
activity begins in the earliest phases of stellar evolution, prior to the arrival on the Main 
Sequence (MS). In many cases, strong and rapid line variability has been detected with 
periods close to the star rotation period, interpreted as due to the rotational modulation of 
spots on the stellar surface or of fast and slow streams in the stellar wind (e.g. Gahm et 
al. 1995; Catala et al. 1989). On the other hand, the fact that young stars during their 
evolution to the MS move across the instability region of post-MS stars raises the possibility 
that at least part of the activity could also be due to stellar pulsations (see Baade & Stahl 
1989; Kurtz & Marang 1995). However, the instability properties of young stars have not 
been studied so far and it is not at all clear whether a Pre-Main-Sequence star can indeed 
pulsate. 

As recently discussed by Gautschy & Saio (1996), stars over essentially the whole 
mass spectrum can become pulsationally unstable during various stages of their evolution. 
In particular, 5 Scuti stars are intermediate-mass variables of spectral type A to F with 
pulsation periods less than 0^3 and light amplitudes ranging from thousandths of magnitude 
to some tenths (Breger 1979; Breger & Pamyatnykh 1998). In spite of the traditional 
association of 5 Scuti pulsation with stars in their core hydrogen-burning phase or evolving 
towards the base of the giant branch and burning hydrogen in a shell, some S Scuti have 
been identified to be PMS objects (Breger 1972, Kurtz & Marang 1995). Kurtz & Marang 
(1995) have observed the Herbig Ae star HR 5999 searching for pulsational photometric 
light variations characteristic of 8 Scuti stars. They have found a peak-to-peak pulsation 
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amplitude of about 013 in Johnson V-band and a period of 4.99 hr, superimposed on a 
nonperiodic photometric variation of m 35 probably due to variations of dust obscuration in 
the circumstellar envelope and disk. There are other stars in which evidence for pulsational 
variability has been advocated to explain the observed rapid line profile variability (e.g. 
HD 163296 and HD 250550, Baade & Stahl 1989), but the case is not as convincing as that 
provided by HR 5999. In any case, the discovery of pulsation in PMS stars is extremely 
important since it provides a unique means for constraining the internal structure of young 
stars and for testing evolutionary models. 

In this Letter, we explore the pulsational behavior of stellar models distributed along 
the evolutionary tracks of PMS stars computed by Palla & Stahler (1993). We identify the 
boundaries of the instability strip by using nonlinear models. We can then use these models 
to constrain the fundamental stellar parameters of HR 5999 and to select a sample of known 
PMS stars that are likely candidates for showing pulsation variability in their photometric 
light curves. 

2. Pre-main-sequence evolutionary models 

The calculation of the pulsational instability requires the specification of the stellar 
parameters, luminosity and effective temperature, for a given mass. We have used the PMS 
models for low- and intermediate-mass stars computed starting at the birthline determined 
by the protostellar accretion phase (Palla & Stahler 1990, 1993). The development of 
protostar theory has altered the concept of PMS evolution particularly for stars more 
massive than solar. In comparison with classical evolutionary tracks that used arbitrary 
initial conditions, those of Palla & Stahler (1993) occupy a much reduced portion of the 
HR diagram as a consequence of the modest stellar radii attained by each star during the 
accretion period. In addition, stars between 2 and 4 M are not initially convective from 
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surface cooling, but begin their evolution in a thermally unrelaxed state corresponding to 
luminosities and effective temperatures lower than those of classical evolutionary tracks. 
This sequence of events may have direct effects on the study of the instability properties of 
PMS stars and should be taken into account in the selection of the equilibrium models. 

As an illustration of the differences in the internal structure of stars at various stages 
of evolution, Fig. 1 shows the variation with radius of the interior density (normalized to 
the central value) for a 3 M Q star. On the birthline, at the beginning of the PMS phase, the 
density profile is quite flat and the star is thermally unrelaxed. At the end of the relaxation 
phase, in the middle of the PMS track, the star has already started its contraction to the 
ZAMS and is more centrally condensed. The middle curve of Fig. 1 is for a model with 
luminosity L*=57 L and effective temperature T e ff=6900 K at an age of 7 x fO 5 yr. As we 
shall see below, such a model falls within the boundaries of the instability strip. Finally, the 
lower curve is for a star that has completed the main-sequence phase and is evolving toward 
the red giant branch. The effective temperature is approximately the same as for the PMS 
star, but the luminosity is now L*=135 L Q due to the larger radius (8.3 R© instead of 5.3 
R©). These conditions correspond to those of a typical, evolved 5 Scuti star. The star is 
much more centrally condensed because of the complete exhaustion of hydrogen. However, 
beyond the inert core, the density profile is quite similar to that of the PMS star. Thus, a 
mature 5 Scuti star differs from a young one mostly in the inner regions and for its higher 
luminosity. These structural changes clearly affect the estimate of the pulsation periods, as 
we shall see in the following sections. 

3. Pulsational models 

The theoretical framework for the investigation of the linear and nonlinear pulsation 
characteristics has been described in several papers (e.g. Bono et al. 1997 and references 
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therein) and is briefly summarized here. We have computed several sequences of linear 
nonadiabatic models at fixed mass covering a wide range of luminosities and effective 
temperatures. These models give information about periods and modal stability and 
provide the static envelope structures for the nonlinear computations. We adopt an optical 
depth of the outermost zone of about 0.001 and fix the inner boundary condition at a 
fraction of the stellar radius varying between 10 and 20%. The zone closest to the Hydrogen 
Ionization Region (HIR) was constrained to Thir = 1.3 x 10 4 K and 20 mass zones were 
inserted between the HIR and the surface in order to ensure a good spatial resolution of the 
outermost regions throughout the pulsation cycle. 

Due to the key role played by spatial resolution for firmly estimating both linear and 
nonlinear modal stability of higher modes, the envelopes were discretized by adopting a 
detailed zoning in mass. The mass ratio between consecutive zones (6) has been assumed 
equal to 0=1.1 for temperatures lower than 6 x 10 5 K and to 6=1.2 for higher temperatures. 
In this way a typical envelope model, used for both linear and nonlinear computations, is 
characterized by roughly 20-30% of the total stellar mass and by 150-250 zones. 

To model the physical mechanisms which govern the pulsation phenomenon, the k 
and 7 mechanisms in the hydrogen and helium ionization regions (see Cox 1980), we 
have adopted the OPAL radiative opacities provided by Rogers & Iglesias (1992) for 
temperatures higher than 10 4 K and the molecular opacities of Alexander & Ferguson 
(1994) for lower temperatures. A solar chemical composition (Y=0.28, Z=0.02) is assumed 
in all the computations. 
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4. Linear analysis 

For each selected pulsation mode, the fundamental observables that can be derived 
from the linear nonadiabatic computations are periods and growth rates. In this work, we 
limited the study to the first three radial modes of pulsation. A detailed analysis of the sign 
of growth rates provides an estimate of the location of the blue boundaries of the instability 
strips for each mode, representing the maximum effective temperatures allowing pulsation 
at a given luminosity level. No information is derived about the location of red boundaries 
because our linear models are completely radiative while the quenching of pulsation at lower 
effective temperatures is due to the interaction between convective and dynamical motions 
(e.g. Baker & Kippenhan 1965). 

Although our goal is the definition of the theoretical instability strip, we have specialized 
the linear analysis to the case of HR 5999. In particular, we begin by reproducing its 
observed pulsation period of 4.99 hr for a combination of mass, luminosity and effective 
temperature consistent with both observations and evolutionary tracks. To this purpose, 
we have computed a grid of linear nonadiabatic models of stars with masses equal to 1.5, 
2, 2.5, 3, 3.5, and 4 M , luminosities in the range 10 to 140 L , and effective temperatures 
4500<T cff <9500 K. The linear pulsation code assumes that the stellar models are in 
complete thermal and hydrostatic equilibrium. 

Figure 2 shows the location in the HR diagram of the lines of constant period (P = 4.99 
hr) for the three pulsation modes and for different stellar masses. At a given mass, each 
line interesects the corresponding evolutionary track in one point that determines the 
values of L* and T e jj to be used in the nonlinear models. Also shown in Fig. 2 is the 
position of HR 5999 with the observational uncertainties. The distance to the star based 
on Hipparcos parallaxes is d = 210t|o P c which gives a luminosity of 78±23 L Q (van den 
Ancker et al. 1997). The spectral classification of HR 5999 is A7III/IVe and estimates of 
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the effective temperature for this spectral type vary between ~7600 K and ~7000 K (e.g. 
van den Ancker et al. 1997; Solano & Fernley 1997). As one can see from Fig. 2, only few 
equilibrium models with the correct period fall within the HR 5999 box. In particular, the 
model with mass M*=3 M pulsates in the fundamental mode, that with M*=3.5 M in 
the first overtone, and that with M*=4 M Q in the second overtone. The selection of these 
models restricts the choice of acceptable stellar masses for the detailed nonlinear analysis. 

5. Nonlinear analysis 

Given the static envelope structures produced by the linear computations, we have 
performed the nonlinear stability analysis. The equations governing both the dynamical 
and convective structures were integrated in time until the initial perturbations and the 
nonlinear fluctuations, which result from the superposition of higher order modes, settled 
down (for more details see Bono, Castellani & Stellingwerf 1995). Because of the very 
low growth rates (between 1CT 5 and 10~ 4 ), an extremely large period number (exceeding 
10,000) must be calculated for each model before reaching the stable limit cycle of pulsation 
(maximum amplitude). The initial velocity profile was obtained by perturbing the linear 
radial eigenf unctions with starting surface velocities varying in the range 0.1 to 10 km s -1 . 

As found in investigations of Post-Main-Sequence 5 Scuti star models (Milligan & 
Carson 1992; Bono et al. 1997), the modal stability analysis reveals the presence of 
stable nonlinear limit cycles in the fundamental, first and second overtone as the effective 
temperature increases. Figure 3 shows the location of the theoretical instability strip 
boundaries in the HR diagram, together with the last static and the first pulsating models 
found for different masses. The stellar parameters along the strip are listed in Table 1. 
In agreement with the distribution originally suggested by Breger & Bregman (1975) and 
whith previous theoretical predictions (Bono et al. 1997), we find that second and first 
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overtone pulsators are located at effective temperatures higher than the fundamental ones. 
All the models on the red edge are pulsating in the fundamental mode, while those on 
the blue edge in either of the overtone modes. The red edge varies between ~6500 K 
and ~7100 K, whereas the blue edge between ~7100 K and ~7500 K. The width of the 
instability strip is nearly constant and equal to about 650 K. The time spent by each star 
within the strip is also a constant fraction of the total PMS contraction time, typically 
between 5 and 10%, although in absolute terms it decreases from ~10 6 yr for M*=1.5 M Q 
to 8xl0 4 yr for M*=4 M . 

As shown in Fig. 3, the PMS instability strip terminates abruptly at the blue edge of 
the M*=4 M evolutionary track. This results is the consequence of our assumption that 
PMS stars begin their contraction phase at the birthline computed with an accretion rate 
of M acc = 10~ 5 M yr -1 . It has been shown that such birthline reproduces quite well the 
the distribution of known PMS stars in the HR diagram (Palla & Stahler 1990). Obviously, 
the theoretical instability strip extends to luminosities higher than the upper limit found 
here for PMS stars, but it would then be populated only by more evolved objects such as 
the classical S Scuti stars. 

6. Discussion 

The models presented in this paper represent the first attempt to derive the boundaries 
of the instability strip for young stars. A complete analysis of the predicted light curves and 
amplitudes and of the sensitivity to the input physics will be given elsewhere (Marconi et 
al. 1998, in preparation). However, we can use these initial results to set constraints on the 
pulsational properties of known PMS objects. The location of HR 5999 with respect to the 
instability strip is shown in panel (b) of Fig. 3. We can see that in order to fall within the 
strip, the effective temperature of HR 5999 should be lower than ~7250 K if its mass is 3 
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M Q or equal to ~7100 K for M*=4 M . However, if we accept the pulsation period of 4.99 
hr found by Kurtz & Marang (1995), then the only acceptable models would be (see Fig. 2) 
a 3 M star pulsating in the fundamental mode, a 3.5 M in the first overtone, or a 4 M 
star in the second overtone. On the other hand, these three models are located at the blue 
edge of the strip and the nonlinear stability analysis predicts that they should all pulsate 
in the second overtone mode. But in this higher mode only the 4 M model oscillates with 
the correct period of HR 5999, ruling out the 3 M and 3.5 M stars as acceptable models. 
Thus, we conclude that the pulsation properties of HR 5999 suggest a mass of 4 M with 
second overtone pulsation. 

Despite the relatively short time spent by PMS stars inside the strip, Fig. 3 shows 
that there is a number of known objects with the right combination of luminosity and 
temperature to qualify as candidate pulsators. These include five stars (V376 Cas, CQ Tau, 
LkHa 208, HD 142666 and KK Oph) of the UX Orionis class, characterized by large 
photometric and polarimetric variations attributed to variable extinction by circumstellar 
dust (Grinin & Tambovtseva 1995; Natta et al. 1997). As in the case of HR 5999, these 
stars should show periodic variations of smaller amplitudes in their photometric behavior. 
Also, there are eight Herbig Ae/Be stars of spectral type between A7 and F7 (HD 35929, 
V351 Ori, HD 142527, LkHa 233, TY CrA, RCr A, TCr A and HD 104237) which are 
located inside or near the instability strip within the uncertainties of the determination of 
L* and T e ff (Berrilli et al. 1992; van den Ancker et al. 1998). Searches for 5 Scuti type 
variability in these stars are now under way. 

It is a pleasure to thank G. Bono and L. Terranegra for useful discussions on pulsation 
models and spectral properties of young stars. This project was partly supported by CNR 
grant 97.00018.CT02 and by ASI grant ARS 96-66 to the Osservatorio di Arcetri. 
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Fig. 1. — The internal density profile (normalized to the central value) of a 3 M star in 
different evolutionary stages. The upper curve is a protostar on the birthline computed with 
M acc =fCT 5 M yr" 1 (Palla & Stahler 1990). The middle curve is for a pre-MS star with 
luminosity L*=55 L , effective temperature T e fj=6900 K, radius R*=5.3 R and an age of 
tpMS = 7 x 10 5 yr. Finally, the lower curve is for a post-MS star with the same effective 
temperature, but L*=135 L and R*=8.3 R . 

Fig. 2. — Location of the lines of constant period (equal to 4.99 hr) in the HR diagram. 
Three sets of lines are shown for the fundamental, first and second overtone of stars of 
different masses. The four curves of the fundamental mode are computed for M*=2.5, 3, 
3.5 and 4 M . The three curves for the first overtone refer to M*=3, 3.5 and 4 M . The 
last two curves for the second overtone are for M*=3. 5 and 4 M . The filled dots mark the 
intersection of the constant period lines with the corresponding evolutionary tracks of Palla 
& Stahler (1993), labelled by the appropriate mass. The long dashed curve is the stellar 
birthline computed with an accretion rate of M acc = 10~ 5 M yr -1 . Finally, the position of 
HR 5999 is indicated by the two asterisks, together with the observational uncertainties in 
L* and T cff . 

Fig. 3. — Left panel: location of the instability strip of PMS stars in the HR diagram. Right 
panel: distribution of known PMS stars within or near the boundaries of the instability strip 
which can be considered candidates for pulsational instability, triangles show the position 
of five UX Orionis-type stars (from Natta et al. 1997); squares are eight Herbig Ae/Be 
stars (from Berrilli et al. 1992; van den Ancker 1998); the pentagon is for HR 5999 and its 
associated uncertainty. 
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Table 1. Parameters of the instability strip 



Red Edge Blue Edge 



Mass 


L* 


T ff 


Mode 


Period 




T ff 


Mode 


Period 


(M ) 


(Lo) 


(K) 




(d) 


(Lo) 


(K) 




(d) 


1.5 


7 


7090 


f 


0.0611 










2.0 


18 


6870 


f 


0.1195 


21 


7500 


II 


0.0690 


2.5 


32 


6640 


f 


0.2155 


37 


7350 


II 


0.0925 


3.0 


52 


6500 


f 


0.2601 


62 


7260 


II 


0.1299 


3.5 


77 


6550 


f 


0.3195 


95 


7200 


II 


0.1721 


4.0 










133 


7100 


I 


0.2736 




0.2 0.4 0.6 0.8 1 



Radius (r/R.) 



